Quantum Storage of a Photonic Polarization Qubit in a Solid 
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We report on the quantum storage and retrieval of photonic polarization quantum bits onto and 
out of a solid state storage device. The qubits are implemented with weak coherent states at the 
single photon level, and are stored for 500 ns in a praseodymium doped crystal with a storage and 
retrieval efficiency of f %, using the atomic frequency comb scheme. We characterize the storage 
by using quantum state tomography, and find that the average conditional fidelity of the retrieved 
qubits exceeds 95 % for a mean photon number fi = 0.4. This is significantly higher than a classical 
benchmark, taking into account the Poissonian statistics and finite memory efficiency, which proves 
that our device functions as a quantum storage device for polarization qubits, even if tested with 
weak coherent states. These results extend the storage capabilities of solid state quantum memories 
to polarization encoding, which is widely used in quantum information science. 

PACS numbers: 03.67.Hk,42.50.Gy,42.50.Md 



The ability to transfer quantum information in a co- 
herent, efficient and reversible way from light to matter 
plays an important role in quantum information science 
[l[ . It enables the realization of photonic quantum memo- 
ries (QM) which are required to render scalable elabo- 
rate quantum protocols involving many probabilistic pro- 
cesses that have to be combined. A prime example is the 
quantum repeater , where quantum information can 
be distributed over very long distances. Other applica- 
tions include quantum networks 6], linear optics quan- 
tum computation Q , deterministic single photon sources 
[U and multiphoton quantum state engineering. 

Proof of principle experiments demonstrating photonic 
QMs have been reported in different atomic systems such 
as cold [9l 4l5l a nd hot atomic gases flB- [l9| . single atoms 
in cavities [20j| and solid state systems |2l[ [22j . In re- 
cent years, solid state atomic ensembles implemented 
with rare-earth doped solids have emerged as a promising 
system to implement QMs. They provide a large num- 
ber of atoms with excellent coherence properties natu- 
rally trapped into a solid state system. In addition, they 
feature a static inhomogeneous broadening that can be 
shaped at will, enabling new storage protocols with en- 
hanced storage properties (e.g. temporal multiplexing) 
23, 24]. Finally, some of the rare-earth doped crystals 
(e.g. praseodymium and europium doped crystals) pos- 
sess ground states with extremely long coherence times 
[25I ] (> seconds), which hold promise for implementing 
long lived solid state quantum memories. 

Recent progress towards solid state QMs include the 
storage of weak coherent pulses at the single photon level 
21, 26-28], the quantum storage of coherent pulses with 
efficiency up to 70 % [22j , the spin state storage of bright 
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coherent pulses [25j, |29j and the storage of multiple tem- 
poral modes in one crystal [3(1 [M| . Very recently, these 
capabilities have been extended to the storage of non- 
classical light generated by spontaneous down conversion, 
leading to the entanglement between one photon and one 
collective atomic excitation stored in the crystal [32I . [H| , 
and entanglement between two crystals [34] ■ 

All previous experiments with solid state QMs have 
been so far limited to the storage of multiple modes us- 
ing the time degree of freedom, e.g. time bin or energy 
time qubits. However, quantum information is very often 
encoded in the polarization states of photons, which pro- 
vide an easy way to manipulate and analyze the qubits. 
Extending the storage capabilities of solid state QMs to 
polarization encoded qubits would thus bring much more 
flexibility to this kind of interface. Unfortunately, stor- 
ing coherently polarization states is not straightforward 
in rare-earth doped crystals. The main difficulty is that 
these crystals have in general a strongly polarization de- 
pendent absorption. Storing directly a polarization qubit 
in such a system would result in a severely degraded fi- 
delity of the retrieved qubits. 

In this paper, we report on the storage and retrieval 
of a photonic polarization qubit into and out of a solid 
state quantum storage device with high conditional fi- 
delity. The photonic qubits are implemented with weak 
coherent pulses of light, with a mean photon number u 
from 0.01 to 36. We measure the conditional fidelity Q 
of the storage and retrieval process (i.e. assuming that a 
photon was re-emitted) and compare it to classical bench- 
marks. With this procedure, we can show that our crys- 
tal behaves as a quantum storage device, even if tested 
with classical, weak coherent pulses. We overcome the 
difficulty of anisotropic absorption by splitting the po- 
larization components of the qubit and storing them in 
two spatially separated ensembles within the same crys- 
tal dimGiil. 

Our memory is implemented using a 3 mm long 
Pr 3+ :Y 2 Si0 5 crystal (0.05%). The relevant atomic 
transition connects the 3 H4 ground state to the 1 D 2 ex- 
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FIG. 1: (Color online) (a) Experimental setup. The prepa- 
ration beam and the qubit beam are derived from the same 
laser at 606 nm. Both beams are amplitude and frequency 
modulated using an acousto-optic modulator (AOM). The po- 
larization of the preparation beam is set to 45 degrees, while 
the polarization of the qubit beam can be set arbitrarily with 
a half wave plate (HWP) and a quarter wave plate (QWP). 
The qubit beam is attenuated down to the single photon level 
using neutral density filters (NDF). The two beams are re- 
combined at a beam splitter (BS) and sent to the storage 
device which is composed of a Pr 3+ : Y2SiOs crystal cooled 
to 2.8 K, two beam displacer (BD) and two HWP. The light 
released by the memory is then sent to a polarization ana- 
lyzer composed of a HWP, a QWP and a polarization beam 
splitter (PBS), before being detected with a single photon de- 
tector (SPD). The two mechanical shutters (SI and S2) are 
used to suppress optical noise from the preparation beam and 
to protect the SPD. (b) Storage and retrieval of a weak V) 
qubit with duration 140 ns and \i — 0.4. The temporal his- 
togram of the detection without (dotted line,empty pit) and 
with (solid line) AFC is shown. The dashed lines around the 
AFC echo show the detection window. Inset: Level scheme of 
Pr 3+ : Y 2 SiO s with the dashed arrows showing the relevant 
transition for photon absorption and re-emission. 



cited state and has a wavelength of 605.977 nm. Each 
state has three hyperfine sublevels as shown in Fig. [T] 
The measured maximal optical depth at the center of the 
5 GHz inhomogeneous line is OD = 7. We use the Atomic 
Frequency Comb (AFC) scheme to store and retrieve the 



qubits 24]. This requires to shape the inhomogeneous 
absorption profile into a series of periodic and narrow 
absorbing peaks, placed in a wide transparency window. 
This creates a frequency grating and when a photon is 
absorbed by the comb, it will be diffracted in time and 
re-emitted after a pre-determined time t$ = 1/A, where 
A is the spacing between absorption peaks. 

Our experimental apparatus is described in Fig. Q] 
The laser source to generate light at 606 nm is based on 
sum frequency generation (SFG) in a PP-KTP waveg- 
uide (AdVR corp) from two amplified laser diodes at 
1570 nm (Toptica, DL 100 and Keopsys fiber amplifier) 
and 987nm(Toptica, TA pro). With input power of 
370 mW and 750 mW for the 987 nm laser and 1570 nm 
laser, respectively, we achieve an output power of 90 mW 
at 606 nm. Taking into account the 30 % coupling effi- 
ciency of both beams into the waveguide, we obtain a 
SFG efficiency of - 350 %W _1 . The beam is then split 



in two parts, one which will be used for the memory 
preparation (preparation beam) and one to prepare the 
weak pulses to be stored (qubit beam). In each path, 
the amplitude of the light is modulated with an acousto- 
optic modulator (AOM) in a double pass configuration, 
in order to create the required sequence of pulses for 
the preparation of the memory and of the polarization 
qubits. The radio-frequency signals used to drive the 
AOMs are generated by an arbitrary waveform generator 
(500Msample/s, 200 MHz, 1GB internal memory, PXIe 
module and ProcessFlow software from Signadyne). Af- 
ter the AOMs, both beams are coupled to a polarization 
maintaining optical fiber and sent to another optical ta- 
ble where the cryostat is located. 

The crystal is cooled down to 2.8 K in a cryofree 
cooler (Oxford Instruments V14). After the fibers, the 
preparation beam is collimated to a diameter of around 
600 fim with a telescope and sent to the storage device. 
Right before the cryostat, a beam displacer (BD1) splits 
the two polarization components of the incoming light 
onto two parallel spatial modes separated by 2.7 mm, co- 
propagating through the crystal. To ensure equal power 
in both spatial modes, the polarization of the preparation 
beam is set to 45 degrees. After BD1, the polarization 
of the horizontal beam (lower beam in Fig[lJ is rotated 
by 90 degrees using a HWP such that both beams enter 
the crystal with the same polarization which is parallel 
to the D2 axis of the crystal, which maximizes the ab- 
sorption. The qubit beam is strongly attenuated by a 
set of fixed and variable neutral density filters (NDF), 
and \i before BD1 was varied from 0.01 to 36. After the 
NDF, arbitrary polarization qubits are prepared, using 
a quarter (QWP) and a half wave plate (HWP). The 
qubits are then overlapped to the preparation beam at 
a beam splitter (BS). After the cryostat, we rotate back 
the polarization of the lower beam and the two spatial 
modes are combined again at a second beam displacer 
(BD2). The two path between BD1 and BD2 form an 
interferometer with very high passive stability [l2l |36| . 

After the interferometer, the transmitted and retrieved 
light enters the polarization analysis stage, composed of 
a QWP, a HWP and a polarization beam splitter (PBS), 
which allow us to measure the polarization in any basis. 
The transmitted beam at the PBS is coupled in a mul- 
timode fiber and sent to a Silicon avalanche photodiodc 
Single Photon detector (SPD, model Count, Laser Com- 
ponents). The electronic signal from the SPD is finally 
sent to a time stamping card (PXIe card from Signa- 
dyne) in order to record the arrival time histogram. The 
mean photon number /j, is determined by measuring the 
detection probability per pulse pdet when no atoms are 
present (i.e. with the laser 60 GHz off resonance), and 
backpropagating before BD1 taking into account the de- 
tection efficiency (tjd — 50 %) and the transmission from 
before BD1 to the detector (rj t = 40%). The prepara- 
tion beam and the qubits are sent sequentially towards 
the crystal. The total experimental sequence lasts 3 sec- 
onds, during which the preparation lasts 1 s. During the 
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FIG. 2: (Color online) (a) Measured detection probability 
(Pdet) as a function of the polarizer angle, for \V) (filled cir- 
cles) and \D) (open circles) input polarization qubits, with 
/i = 0.4. The fitted raw visibilities are (97 ± 0.5) % and 
(83 ± 2) %, respectively, (b) pdet as a function of polarization 
angle for \R) polarization qubit input, with (filled squares, 
(89 ± 2) %) and without (open squares, (15 ± 3) %) QWP in- 
serted before the polarizer. 
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FIG. 3: (Color online) Quantum state tomography. Recon- 
structed density matrices of the retrieved qubits for \H), \D) 
and |_R) input qubits, with /i = 0.4 . No background has been 
subtracted. 



next 2 seconds 10 5 weak pulses are prepared, stored and 
retrieved at a rate of 50 kHz. 

In order to create the AFC, we first create a wide trans- 
parency window within the 5 GHz inhomogeneous profile. 
This is achieved by sending a series of pulses of duration 
1.1ms, during which the frequency of the light is swept 
linearly over a range of 12 MHz. The AFC is then cre- 
ated using the burn back procedure introduced in Ref. 
[37l ], i.e. by illuminating the sample with short pulses of 
duration 2 ms, while shifting the frequency of the light by 
—27 MHz with respect to the center of the pit. Four burn 
back pulses are sent with different frequencies separated 
by the comb spacing, leading to a 4-tooth comb. 

As a first experiment, we verified that a complete set 
of qubit distributed over the Poincare sphere could be 
stored and retrieved in the AFC. We set the storage 



Input State Fidelity Input State Fidelity 
\H) 0.982 ± 0.003 \V) 0.983 ± 0.002 
\D) 0.968 ± 0.005 \A) 0.938 ± 0.009 
\R) 0.954 ± 0.007 \L) 0.926 ± 0.01 



TABLE I: Raw conditional fidelities for 6 different polariza- 
tion input states. For this measurement, T}m varied between 
8 % and 10 %. The errors have been obtained using Monte 
Carlo simulation taking into account the statistical uncer- 
tainty of photon counts and a technical error reflecting slow 
drifts in our systems, and estimated from the residuals from 
the fit of Fig[2] and similar curves. 



time to 500 ns and we used the following input states: 
\H),\V), \D) = ±(\H) + \A) = -±(\H) \V)), 

\ R ) = + ^ and l L > = 75<-W i\V))- FiglD 

(b) shows the experimental storage of a \V) qubit en- 
coded onto a pulse of duration 140 ns (FWHM) and 
with fi = 0.4. Similar curves are obtained for the other 
states, and the average storage and retrieval efficiency is 
T)M — (10.6±2.3)%. In order to test that the coherence 
between the \H) and \ V) components of the qubits is pre- 
served during the storage and retrieval, we then recorded 
the number of counts in the retrieved pulses when rotat- 
ing the detection polarization basis using a HWP, for var- 
ious input states. In Fig. H](a), the curves obtained for 
\V) and \D) are shown. We observe interference fringes, 
with a raw fitted visibility of (97 ± 0.5)% for the \V) 
qubit, and (83 ± 2) % for the \D) qubit. In the case of a 
perfect circular \R) qubit, we should observe no depen- 
dance on the HWP angle when no QWP is inserted. The 
interference should be restored however, when a QWP 
is inserted before the HWP, which turns circular polar- 
ization into a linear one. The results are shown in Fig. 
[H (b). We indeed observe a strongly reduced visibility 
without the HWP ((15 ± 3)%), while a fringe with a 
high visibility of (89 ± 2) % is obtained with the QWP. 
The residual visibility without the QWP may be due to 
a non perfect preparation of the \R) state. The non per- 
fect visibility for the \D) and \R) is mostly due to small 
phase fluctuations engendered by mechanical vibrations 
from the cryostat. Indeed, we observe similar visibilities 
for bright pulses out of resonance with the atomic transi- 
tion. These results show that the phase between the two 
polarizations components is almost perfectly preserved in 
the storage and retrieval process, and is preserved to a 
high degree in our combined interferometer and memory 
setup, for various polarization input states. 

In order to better characterize the quality of the stor- 
age process, we reconstruct the density matrix of the 
retrieved qubits using quantum state tomography [38| . 
for the complete set of qubits described above, with 
fx = 0.4. The reconstructed output density matrices 
Pout for \H), \D) and \R) are shown in Fig [3J From 
the matrices p ou t, we can then estimate the conditional 
fidelity of the output states with respect to the target 
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FIG. 4: (Color online) Average fidelity measured as a func- 
tion of the mean number of photon per pulse fj,. The fidelity is 
measured by quantum state tomography and is the average of 
3 input states V), \D), and \R). Filled circles are experimen- 
tal data without any background subtraction. Empty squares 
correspond to dark count subtracted fidelities. The error takes 
into account statistical uncertainty of photon counts, techni- 
cal errors and standard deviation of fidelities for the 3 polar- 
izations. The various lines correspond to classical thresholds 
for different situations. The horizontal line is the limit of 2/3 
for single qubits (N = l). The dashed-dotted line corresponds 
to Eq[T] where the Poissonian distribution of photon number 
is taken into account. Finally, the two other lines correspond 
to the cases where the finite storage efficiency is taken into ac- 
count. The solid line corresponds to r\ — (10±2) %. Measured 
r\M are between 8 and 10 % for all points. The dashed line 
corresponds to r\ — r]MT]tT]D = 2%. For fi > 1, an additional 
ND filter is placed before the detector to avoid saturation 
effects. 



state .Fj^A = (ip\pout W)- The values for the complete set 
of inputs are listed in Table 1. We find a mean fidelity 
of F^ ean = (96 ± 2) %. We emphasize that this value is 
a lower bound for the conditional fidelity, since it is cal- 
culated with respect to a target state and also takes into 
account imperfections in the preparation of the qubits. 

Finally, in order to assess the quantum nature of the 
storage, we determine the average fidelity as a function 
of n, and compare it with the best obtainable fidelity 
using a purely classical method consisting of measuring 
the state and storing the result in a classical memory. It 
has been shown that for a state containing N qubits, the 
best classical strategy leads to a fidelity of F c — (N + 
l)/(N + 2) HU, leading to the well known fidelity of 2/3 
for N = 1. If the qubit is encoded in a weak coherent 
state, as it is the case in our experiment, one has to take 
into account the Poissonian statistics of the number of 
photons [13], and the classical fidelity is given by: 



F 



class 



JV + l P(p,N) 



N>1 



N + 2 1-P(^,0) : 



(1) 



where -P(/i, N) = e~^{i N /N\. This is valid for the case 
of a memory with unity efficiency. If t}m < 1 > the classi- 
cal memory could use a more elaborate strategy to take 
advantage of finite efficiency in order to gain more infor- 



mation about the input quantum state [2Cf . For example 
the classical memory could give an output only when the 
number of photons per pulse is high, and hence estimate 
with better fidelity the quantum state (See appendix [A"|. 
The different curves corresponding to the discussed cases 
are plotted in Fig. |4] as a function of fi. The points 
correspond to experimental data. Measured F^ ean are 
significantly higher than the classical fidelity, for most of 
the photon numbers tested. This proves that our device 
performs as a quantum storage device for polarization 
qubits, even if tested with weak coherent states. We ob- 
serve that the measured raw fidelity decreases for /i < 0.1. 
This is mainly due to the dark count of the SPD, as high 
fidelities can be recovered by subtraction of this back- 
ground (open squares). We also observe that when /i 
becomes too large (fi > 3.5 in our case), the measured fi- 
delity is not sufficient to be in the quantum regime. This 
confirms that very low photon numbers are required to 
test the quantum character of QMs with weak coherent 
states. To our knowledge, it is the first time that an en- 
semble based memory has been characterized using this 
criteria. 

The storage time in our experiment is limited by the 
minimal achievable width of the AFC peaks (<600 kHz), 
which is in turn fully limited by the linewidth of our un- 
stabilized laser. Peaks as narrow as 30 kHz have been cre- 
ated in Pr 3+ : Y 2 Si05 using a frequency stabilized laser 
J40|, which should allow a storage time in the excited 
state of about 10 /is. This should also allow the storage 
of multiple polarization qubits in the time domain. In or- 
der to increase the storage time and to achieve on demand 
read out with an AFC, the optical excitations should be 
transferred to long lived spin excitation as demonstrated 
for bright pulses in 29]. 

We have demonstrated the quantum storage and re- 
trieval of polarization qubits implemented with weak co- 
herent pulses at the single photon level, in a solid state 
storage device. The conditional fidelity of the storage 
and retrieval is > 95 %, significantly exceeding the clas- 
sical benchmark calculated for weak coherent pulses and 
finite memory efficiency. We thus show that solid state 
QMs are compatible with photonic polarization qubits, 
which are widely used in quantum information science. 
This significantly extends the storage capabilities of these 
types of memories. By combining the time and polar- 
ization degrees of freedom one could readily double the 
number of modes that can be stored in the memory and 
create quantum registers for polarization qubits. Using 
these resources, it may also be possible to design a quan- 
tum memory for complex light states such as hyperentan- 
gled states. 

We note that related results have been obtained by two 
other groups 
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Appendix A: Conditional Fidelity using Classical 
State Estimation for Weak Coherent States 



In this section, we give more details on the calcula- 
tion of the classical benchmark presented in Fig. 4 of 
the main manuscript. The general idea is to estimate 
what is the best efficiency that can be obtained using a 
classical method. In particular, we consider a measure 
and prepare strategy where the user performs a classical 
state estimation on the input qubit, stores the result in 
a classical memory and prepares a new qubit according 
to the result obtained. The maximum achievable classi- 
cal fidelity for a state with a fixed photon number N is 
known to be 13911 



F 



N + l 
N + 2' 



If one tests a memory with a true single photon input 
then the classical bound is 2/3. If one uses weak coher- 
ent states then one has to consider the finite probability 
of having more than one photon. A pulse of light with 
mean photon number \x has a probability distribution of 
a Poissonian, given by 

u N 

P( fX ,N) = e-^. 

Then, the maximum achievable fidelity becomes a 
weighted sum over N of the fidelity for a given N where 
the weight is given by the Poissonian statistics of the 
input [20(. We state this below as 



F 



E 

JV>1 



N + l P(n,N) 



N + 2 l-P(/i,0) 



(Al) 



This formula is valid if one assumes that the quantum 
memory has a storage and retrieval efficiency tjm = 1- 
If rjjvf < 1, a more sophisticated classical strategy could 
simulate non-unit efficiency by only giving a result for 
high photon number N, resulting in a higher achievable 
fidelity, as suggested in [2Cj. We consider that the clas- 
sical measure and prepare strategy has an efficiency of 
1, but we define the effective classical efficiency r/c as 
the probability that the classical device gives an output 
qubit, if it has received at least one photon as input: 



Vc 



i-p(m,o) 



where 



Pout — 



E 

N>N min +l 



P(fi,N). 



Note that the photon number statistics of the output 
qubit is not relevant in our case, since we use non photon 



number resolving detectors. The classical memory gives 
a result for some threshold photon number N m - ln + 1 and 
above, and no result for N lower than this. It is impor- 
tant to note that for the above there exists only certain 
r]c for a given mean photon number /i. Hence, not all 
quantum efficiencies can be simulated. A more general 
form of P,ut which allows for arbitrary number and hence 
efficiency is the following 



7 



E p (^ N )- 



AT>iV m i n + l 



where now the memory gives a result for A m j n with a 
probability of 7, with the condition 7 < P(fi, N m i n ). 
The efficiency is then 



VC 



7 +E N >N mia+ i p(m,aQ 
i-p(m,o) 



(A2) 



We now assume that tjm = Vc and N m - m is obtained as 
follows 

A min = mini: PQi,N)<0--P(j*M*lM- (A3) 

N>i+1 

Figure [5] shows a graphical representation of obtaining 
Nmi n . The maximum achievable classical fidelity using 
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FIG. 5: An example of obtaining A m i n . Plotted above is 
the probability distribution for a coherent state of = 3. 
The target efficiency here is 50%. For this target efficiency 
Amin = 3, found using Equation IA3I The grey shaded part 
of bar N = 3 plus the black bars N > 4 equate to the 50% 
target efficiency. 



the above described strategy is then 



AUn + 1'\ \ - 

7+ X. 



N n 



N+l 

N + 2 



P(/i,iV) 



class 



7 + Yl p (^) 

N>N mSn +l 



where 7 and A m j n are obtained from Equations (|A2| IA3|) . 
Figure |5] shows the fidelity as a function of mean photon 
number for various efficiencies 77. 



(A4) 



6 



l 




0.65 



10" 4 10 3 10 2 10 1 10° 10 1 10 2 

FIG. 6: Maximum achievable fidelity using a classical mem- 
ory as a function of mean photon number /i for different 
memory efficiencies 77 (Eq. IA2|) . The cases shown are 
r) = {0.001,0.01,0.1,0.25,0.5,1}. The solid lines are calcu- 
lated using Equation I A4I the dots are calculated using Equa- 
tion [XT] It is seen that for r) = 1, Equation IA4I reduces to 
Equation IA1I The dashed line is for the single photon case 
N = 1. 



Note that the classical memory could also in prin- 
ciple take advantage of the optical loss and the finite 
detection efficiencies in the experiment to increase the 
maximal classical fidelity. In that case, we would have 
Vc = VMVtVD, where r\t is the optical transmission from 
the quantum memory to the detector and t]d is the detec- 
tion efficiency of the SPD. For our experiment, we have 
T]M = 0-1? Vt = 0.4 and rjo — 0.5 such that in that case 
r/c = 0.02. 
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